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Abstract. The formation of the electron energy distribution function in nanosecond 
atmospheric nitrogen discharges is investigated by means of self-consistent solution 
of the chemical kinetics and the Boltzmann equation for free electrons. The post¬ 
discharge phase is followed to few microseconds. The model is formulated in order 
to investigate the role of the cross section set, focusing on the vibrational-excitation 
by electron-impact through resonant channel. Four different cross section sets are 
considered, one based on internally consistent vibrational-excitation calculations which 
extend to the whole vibrational ladder, and the others obtained by applying commonly 
used scaling-laws. 


1. Introduction 

A lot of work has been presented over the last 40 years aimed at understanding the non- 
equilibrium vibrational kinetics of non-equilibrium molecular plasmas Kin 13] due to 
their applications in many technological helds including plasma-assisted combustion jl] 
E], negative ion sources P K E] E]; excimer lasers [10] , microelectronics m, 
aerospace [inila] and environment applications |14] . 

In particular, nitrogen containing plasmas have been extensively considered as 
N 2 is the main component of air. Particular attention has been devoted to the 
coupling of vibrational and electronic distributions with the electron energy distribution 
function through the so-called superelastic collisions (also known as second kind 


2 


collisions) [T51 [TB] . The treatment of non-equilibrium vibrational and electronic kinetics 
has become progressively more complex, with many elementary processes included to 
better model the plasma [171 UHl HSl EHl EH EH EH Elj- Sophisticated models are now 
available which can be used to understand the properties of nitrogen containing plasmas. 
Unfortunately the available experimental measurements generally characterize average 
quantities, such as the degree of dissociation and ionization, rather than the more 
informative vibrational distribution functions (vdf) and electron energy distribution 
functions (eedf), which remain difficult to measure. The recent measurements of the 
vdf of nitrogen [25l [26] are limited to vibrational levels with n < 20 and thus fail to 
probe the plateau predicted for the higher lying vibrational levels. Similarly, the scant 
experimental information on eedf emphasizes the post-discharge regime. Different sets 
of cross sections as well as different forms of non-equilibrium vibrational distributions 
were used by Dyatko et al. [ETIEH] to reproduce experimental vibrational temperatures 
in low pressure nitrogen afterglows. More recently Dyatko et al. have developed a 
complete plasma kinetic model determining the cross sections for electronic excitation 
for n > 0 shifting the known cross sections for the ground vibrational level and making 
use of Franck Condon factors m- A number of theoretical calculations of the vdf and 
eedf are available which are generally characterized by the use of different data sets 
for the key cross sections. In this context interesting results were reported for the low 
pressure nitrogen afterglow by Guerra et al. by coupling eedf and electronically excited 
state kinetics to reproduce the experimental eedf from second derivatives of digitized 
probe characteristics measured using a triple probe technique 1301 . 

The accuracy and the completeness of the existing data bases is not well-known 
and hard to estimate. A sensitivity analysis would be welcome to elucidate the key 
processes which determine the form of both the vdf and eedf. One process is of 
paramount importance for coupling the two distributions, namely the electron-vibration 
(e-V) processes i.e. the pumping of vibrational quanta over the vibrational ladder 
through resonant electron-molecule interaction. For this process a complete set of cross 
sections interconnecting the whole N 2 vibrational ladder through single- and multi¬ 
quantum transitions has recently become available [311 132 ] • These e-V cross sections are 
indeed essential for determining the initial distribution of vibrational quanta over the 
vibrational ladder and also, at sufficiently large electron densities, for giving the spread 
of these quanta over the whole vibrational ladder. In the high electron density case, 
the e-V process plays the same role of the vibration-vibration (V-V) energy transfer 
processes. 

In the past, the lack of a complete cross section database forced researchers to either 
limit their treatment of e-V processes to those interconnecting the hrst 8-10 vibrational 
levels [33l [M] or to use scaling-laws of unknown accuracy to interconnect levels with 
V > 10. Recently, a cross section data set due to Huo et al. [35] was combined with 
extended experimental results to derive a new scaling-law [36] • R should be noted that 
a similar scaling-law was proposed by Gordiets et al. limiting however its application to 
n = 10 [37] • However, all previous models suffer from the incompleteness of the relevant 
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e-V cross section set; this can now be addressed by use of the recent, comprehensive 
compilation of Laporta et al. [311 [32]. 

The e-V processes compete with the V-V and vibration-translation (V-T) energy 
transfer processes in the formation of the non-equilibrium vdf of nitrogen, where the 
V-T process includes as deactivating partners both molecules (V-Tm) and atoms (V- 
Ta). V-V and V-Tm rates have been and are being continuously updated by different 
groups. The htted rates based on the pioneering work by Billing et al. [38l [39] are still 
considered as benchmark values for the numerous other calculations appearing in the 
literature. We note that the hrst studies on the non-equilibrium vibrational kinetics 
including e-V over the hrst few levels, V-V and V-Tm over the whole vibrational ladder 
yielded the well-known Treanor-Gordiets distribution i.e. a vibrational distribution 
containing essentially a Treanor law over the hrst few vibrational levels, followed by 
a long plateau for higher v which hnally ends in a Boltzmann distribution at the gas 
temperature. The interplay between V-V and V-Tm rates determined the end of the 
plateau and the onset of a vdf at the gas temperature Tg. Under favourable conditions 
the long plateau speed-up the dissociation process as well as the ionization one, as 
a result of heavy particle collisions between vibrationally excited molecules (the so- 
called ‘pure vibrational mechanisms’). In these models, V-Ta were constrained to have 
the same rates as V-Tm i.e. were inefficient at destroying the vibrational content of 
nitrogen molecules. This assumption persisted until Lagana et al. [IQI SI] presented 
quasi-classical calculations of the V-Ta processes based on a semi-empirical potential 
energy surface (PES). These data were hrst used by Armenise et al. [12] for describing 
the non-equilibrium vdf of nitrogen under discharge and aerospace re-entry conditions. 
The use of the new rates for discharge conditions yielded to a decrease of the length of the 
plateau in the vibrational distribution due to the interplay of V-V and V-Ta processes. 
The pure vibrational mechanisms were seen to lose their importance in the dissociation 
and ionization processes. A problem however arises when using the V-Ta rates of Lagana 
et al. at low temperature, because these rates were limited to n < 10. This problem has 
been resolved by the extensive quasi-classical trajectory (QCT) calculations on the same 
process by Esposito et al. [IS] [HJ [15] who used the same PES as Lagana et al. These new 
data, while basically conhrming the V-Ta rates of Lagana et al. for v > 10, are orders of 
magnitude lower than the extrapolated Lagana et al. rates for v < 10. These differences 
are extremely important in determining the behaviour of low-lying vibrationally excited 
levels of N 2 which has consequences for the kinetics of the whole system. The use of 
the new e-V and V-Ta cross sections and rates is therefore important to illuminate the 
role of the corresponding processes in affecting the non-equilibrium vibrational kinetics 
of nitrogen plasmas. In particular, they are likely to key for determining structures in 
both the vdf and eedf, with consequences for the dissociation and ionization rates. 

To avoid possible compensation effects with other important elementary processes, 
we base our study on a kinetic code with all other input data hxed at that already 
adopted by our group. Furthermore, to avoid the need for full dimensional fluid 
dynamics, we use a zero-dimensional code based on the coupling of vdf, eedf and 
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electronically excited state kinetics. In particular, we restrict our analysis to the role of 
e-V processes in affecting vdf and eedf of nitrogen discharges i. e. using for the VTa the 
recent analytical curves of Esposito et al. [43] . 

The aim of this paper is to apply these ideas to nanosecond discharges followed 
by a post-discharge analysis of about 1 /rs. In particular, we focus our attention 
on the formation of the eedf under discharge and post-discharge conditions, and the 
creation of structures in the eedf due to superelastic collisions from vibrationally and 
electronically excited states. Particular emphasis is given to the dependence of computed 
microscopic and macroscopic quantities on the electron-molecule resonant vibrational- 
excitation cross sections. The paper is divided in three sections: The Section [2] deals 
with a comparison of the complete e-V data base with approximate scaling laws given 
in the literature. The Section |3] is computation of the eedf and vdf in the nanosecond 
and microsecond regimes, and its dependence on the e-V cross sections. Conclusions 
and perspectives are reported in the Section 01 

2. Theoretical electron-nitrogen cross sections 

Laporta et al. (an [32] calculated low-energy resonant cross sections and corresponding 
rate constants for electron-N 2 . Vibrational-excitation transitions between all vibrational 
levels of ground state N 2 , parameterized on rotational quantum number J, with AJ = 0, 
were considered, using the following scheme: 

e + N 2 (X V, J) ^ N-(2n,) ^ e + N 2 (X ^2+ v', J ), (1) 

where the N 2 (^Ifg) is the short-lived, anionic resonance state. In this section some 
results are reviewed and, in particular, the linear scaling-law between rate constants are 
critically discussed in the light of this new data. 

Laporta et al. used a phenomenological potential energy curve for the N 2 electronic 
ground state which supports 68 vibrational levels for the J = 0 rotation state. The cross 
section calculations were performed using a projection operator formalism within the 
local approximation |16] to solve the nuclear dynamics. Figure [T] shows some results of 
cross sections for single-quantum transitions (upper frame) and for n-quantum transition 
from the n = 0 (middle frame) and from the v = 20 states (lower frame). Typical 
resonant cross sections are characterized by a series of narrow peaks that correspond 
to the vibrational levels of the resonant state. The corresponding rate constants, see 
right-side of Fig. [H are calculated as a function of the electron temperature by the 
convolution of the relevant cross section with a Maxwell distribution for the electrons. 
The complete set of cross sections and rate constants are available from Phys4Entry 
database m- 

As already noted, a complete, internally-consistent, data set of cross sections for 
electron impact vibrational-excitation was not previously available. To complete the 
available data different scaling-laws have been applied. The first scaling-law we consider 
was reported in Ref. [36| and was obtained by fitting the data calculated by Huo et 
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Figure 1. Some series of cross sections (on the left) and the corresponding rate 
coefficients (on the right) for resonant vibrational-excitation of N 2 by electron-impact 
calculated by Laporta et al. isuEa. 


al. [35]. The resulting expression gives rate coefficients for transitions involving excited 
vibrational levels, starting from the rate for the n = 0 level; 
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where a = 0.15. Figure [2] shows the ratio as a function of the vibrational 

level V for different n-quantum jumps at a fixed electron temperature Tg. We note that 
the scaling relation of Eq. (]2|) is a good approximation at high electron temperature, 
but loses accuracy for high vibrational states. 

Another simple scaling-law, reported in Ref. [I9], is given by: 

(*^) ^0,n(^ ) ( 2 ) 

where Ae is the energy difference due to the anharmonicity of the vibrational ladder: 


Ae — Cn — eo + €n+v — Gy . (4) 

Only transitions with v -h n < 10 have been considered. A similar approach has been 
used in Refs. [lailH], where Ae were neglected. 
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Vibrational level v Vibrational level v 

Figure 2. Comparison between the ratio Ky^y+n /and the scaling-law in Eq. (I2|) 
for Tg = 5000 K and Tg = 10000 K, for different vibrational n-quantum jumps, as a 
function of the initial vibrational state. The solid line is the scaling factor 1/(1 -baw). 


A third approach was obtained by extending the scaling-law for V-Tni rates for 
harmonic oscillator to the electron impact rates i.e. it was applied a Laundau-Teller 
model to the e-V rates: 

Kv,v+i = (n + l)A’o,i. (5) 

In the next section we compare results obtained from four different kinetic models which 
use, respectively, the scaling laws in Eqs. (El)-(El) given above and the full set of calculated 
cross sections by Laporta et al. 

3. Kinetic model results 

The kinetic model used in this paper has been described in detail elsewhere [l9l [50l 
[511 |52l [53]. Briefly it consists of homogeneous kinetic equations coupling the non¬ 
equilibrium vibrational and electronically-excited state kinetics with the Boltzmann 
equation for the eedf. The processes induced by heavy particle collisions are listed in 
Table 1 of Ref. [51]. In particular, the processes coupling the vdf with electronically 
excited states and the chemical processes (dissociation and ionization) involving 
electronically excited states should be noted. The list of electron-impact transitions 
with vibrationally excited N 2 molecules is reported in Table 2 of Ref. [5T]. 

The coupling between the chemical kinetics and the Boltzmann equation for 
free electrons is self-consistent, i.e. at each time step the level population and the 
gas composition give the input quantities for the Boltzmann equation and the rate 
coefficients for the electron induced processes are calculated from the actual eedf. 
Inelastic and superelastic collisions involving vibrationally and electronically excited 
states and free electrons are fully considered in this model. Also electron-electron 
collisions are included using an efficient algorithm described in Ref. [5l]. The N 2 
vibrational ladder includes 68 levels. This requires a rescaling of the old V-V and 
V-Tm rates which are based on 46 vibrational levels. 

To solve the coupled problem we impose initial conditions for all the relevant 
quantities. In particular a Boltzmann distribution function at T = 500 K was selected 
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Case Set of cross sections considered 
A Scaling-law of Eq. © is applied to the cross sections [3B] 

B Scaling-law of Eq. ([3]) is used to calculate cross sections up to u = 10 [19] 

C Scaling-law of Eq. m is applied to the cross sections 

D Eull set of cross sections taken from Ref. [SI [32] 

Table 1. Scheme of the models used in the calculations. 

for vdf and a Maxwell distribution function at the same temperature was selected for 
eedf. An initial molar fraction of electrons of Xe = 10“^° was selected while pressure 
and gas temperature were kept at p = 1 bar and Tg = 500 K. Moreover, an uniform 
reduced electric held of E/N = 200 Td was applied for 3 ns and the simulation was 
followed in the post-discharge {E/N = 0) for up to few microseconds i.e. we are using 
an Heavyside-tipe temporal evolution for E/N. The simulation is therefore typical for 
a nanosecond high voltage-high pressure discharge followed by a post-discharge regime 
typically met in plasma-assisted combustion. This kind of simulation can be justihed 
in the framework of the aim of the present paper mainly devoted to the understanding 
of the role of different scaling laws in the e-V cross sections. More realistic E/N forms 
derived from the experiments [5] ES] EH] including a rise time of the applied voltage 
followed by a plateau as well as more complicated forms are left to future work. 

We consider four test cases, see Table [H based on the use of different scaling- 
laws for generating the electron-collision vibrational-excitation cross sections and on the 
numerically calculated data of Laporta et al. Cross sections for all the other process, 
including the heavy particle kinetics, are the same for all the four cases. It should 
be noted that only models A and D include all the transitions between vibrational 
states, while model C considers transitions for levels with n < 10 and single-quantum 
excitations only. 

Before examining the results we estimate the characteristic relaxation times for 
the eedf and vdf. The eedf, neglecting vibrationally excited molecules, reaches a quasi¬ 
stationary state in a timescale of the order, 

tqse = {NKev)~^ ~ 0.3 US , (6) 

where N = 1.52 x 10^® cm“^ is the number density of heavy particles at p = 1 bar and 
T = 500 K and K^v ~ 2 x 10“® cm“^ is the rate coefficient for pumping vibrational 
energy in the system at E/N = 200 Td including all transitions from n = 0. 

The eedf starts to be affected by the presence of vibrationally excited states, as 
well as electronically excited states, when their population starts to be important. 
As an example, the vibrational energy supplied by electron impact is proportional to 
the electron density times the e-V rate; therefore, by considering the maximum value 
obtained for the electron density (ug ~ 1.5 x 10^® cm“^), the characteristic time for 
vibrational excitation is given by: 

-1 


tqsv (^e A"eV ) 




3.3 ns. 


( 7 ) 
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Figure 3. Time evolution of (a) mean electron energy; (b) nitrogen vibrational 
temperature; (c) free electrons molar fraction; (d) nitrogen atoms molar fraction; (e) 
and (f) molar fraction of N 2 triplet states for the four models considered. 


We approximately recover these times in Fig. [3]^a) where the average energy of electrons 
is shown as a function of time for models A to D. 

We start by discussing macroscopic quantities which can be useful for understanding 
the behaviour of the vdf and eedf under discharge and post-discharge conditions. Figure 
Ela), as already noted, reports the time evolution of the electron average energy 
calculated from the eedf as a function of the different models. We can see that in 
the range 0.1-1 ns the mean electron energy remains constant which is typical of an 
eedf in the cold gas approximation, i.e. without the action of superelastic collisions. 
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while the energy undergoes a sudden increase in the time interval 1-3 ns due to the 
excitation of vibrational and electronic states. For t = 3 ns, the electron mean energy 
abruptly decreases as the applied field is switched off. For f > 3 ns, the mean electron 
energy is sustained by the vibrational and electronic energy content of molecules, which 
through superelastic collisions, heat cold electrons. Our results depend strongly on the 
model used for the electron-molecule cross sections, models A and D predict values of the 
average energies lower than the corresponding ones from models B and C. The differences 
in the average energies for the four models can be recovered from the corresponding 
differences in the vibrational temperature calculated from the first two vibrational levels 
(see Fig. |3](b)). We note that case C deviates from the results of cases A, B and D 
mainly due to the poor representation of the vibrational excitation with the scaling- 
law given by Eq. (j5]) which generates much pumped vibrational distribution functions 
(see Fig. Ilj). These differences are important under post-discharge condition but remain 
negligible during the discharge. Inspection of the results shows similarities between 
the mean electron energy obtained from the four models. The vibrational temperature 
initially shows a quasi-stationary value, then grows for nearly 3 ns, where the mean 
electron energy shows a peak. In the post-discharge, the vibrational temperature still 
increases, unlike the mean electron energy, as it is sustained by chemical processes 
(recombination) and electronically excited states. Note that the temporal trend of 1-0 
vibrational temperature is practically the same as the energy stored in the vibrational 
mode. In the long time range (10“'^ — 10“^ s) model B and C show a decrease of the 
vibrational energy, while for models A and D the cooling phase is moved to longer times. 
At longer times the vibrational temperature will converge to the gas temperature. 

Fig. EUc) and (d) report the time evolution of free electron and atomic nitrogen 
molar fractions. Again models A and D give very similar results while models B and 
C predict a much higher concentrations of electrons and nitrogen atoms due to the 
corresponding differences in eedf as well as in the the electron density. In particular 
the electron density is affected by electron impact ionization collisions involving the 
plateau of the vibrational distribution functions which are much more pumped in the 
cases B and C (see Fig. 0]). The lower recombination rates of cases B and C compared 
with the cases A and D are probably due to additional ionization channels promoted by 
the plateau of the vibrational distribution function as well by electron impact collisions 
with electronically excited states. The behaviour of the two triplets of nitrogen molecules 
shown in Fig. [3](e) and (f) is also interesting: their concentration follows the same trend 
as the electron density and reaches a maximum value at the end of the ns discharge. 
Large differences are observed for t > 1 ps where models B and C predict much higher 
concentrations due to more effective recombination. The behaviour of the two triplets 
of nitrogen molecules is also interesting: their concentration follows the same trend as 
the electron density under discharge conditions reaching a maximum value at the end of 
the ns pulse. Large differences are observed for f > 1 ps where models B and C predict 
much higher concentrations as a consequence of 3-body N atom recombination forming 
N 2 (i?) (and then N 2 (A) through radiative decay). Therefore N 2 (A) and N 2 (i?) states 
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Figure 4. Vibrational distributions function of N 2 molecules at three characteristic 
times and for the four models considered in the text. 

follow the behaviour of nitrogen atoms as reported in Fig. |3t^d). 

Figure 0] shows the behaviour of the vdf obtained by the four models at different 
characteristic times: (a) at the end of the discharge pulse; (b) at t = 10 ns; and (c) 
at t = 1 /is in the post-discharge. Models A and D give similar vdfs up to 6 eV, for 
t = 3 and 10 ns, but predict large differences above 6 eV. This behaviour reflects the 
limitations of the scaling-law of Eq. ([2]), as can be observed in Fig. |2l Models B and C 
behave similarly up to 10 ns. At 1 ps, model C predicts a quasi-Boltzmann distribution 
but model B shows a well-developed plateau, in line with models A and D, but at much 
higher v than in these cases. 

Let us now consider the time evolution of eedf under discharge and post discharge 
conditions. Figure [5] illustrates the situation during the discharge. The evolution of the 
eedf towards a quasi-stationary state is rather regular. After 0.1 ns, the eedf reaches 
a quasi-stationary value and, in this regime, superelastic collisions do not affect the 
distribution due to the high value of the electron energy. As shown in Fig. [6l the 
four models give small differences in the eedf, especially in the low-energy region, being 
dominated by e-V inelastic processes. The situation completely changes in the post¬ 
discharge regime where the decrease of average energy enhances the role of superelastic 
collisions from electronic excited states, i.e.: 


e(e) + N 2 (A ^ e(e -f e*) + N 2 (X ^S+, v ), 

(8) 

e(e) + N 2 (B ^n) ^ e(e + e*) + N 2 (X ^S+, v ), 

(9) 

e{e) + N 2 (C ^n) ^ e(e + e*) + N 2 (X v ), 

(10) 


where e* is the threshold energy of the corresponding transition, that are, relative to the 
ground state respectively 6.17 eV, 7.35 eV and 11.03 eV. At f = 10 ns these processes 
are responsible of the peak at e ~ 6 eV, while the second peak at e ~ 12 eV is due 
to electrons in the peak at e ~ 6 eV subjected to another superelastic collision. It 
should be noted that the peaks are quite large due to the superposition of the different 
superelastic collisions. The effect of transition from the C ^11 state are visible as small 
shoulders around 10 eV and 20 eV. These maxima are strongly smoothed by the model 
C due to the higher mean electron energy, which hides the role of superelastic collisions. 
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Figure 5. (a) Time evolution of eedf under discharge conditions according to model 
D. (b) The same plot in (a) in the energy range 0-10 eV. 




Figure 6. (a) eedf at two times during the discharge; (b) full range eedf during the 
post-discharge for the four models. 


Similar results are observed for t = 1 /rs, where only the contribution of the state 
is clearly evident. 

These results clearly demonstrate that the use of complete sets of cross sections, 
as well as the coupling between vibrationally and electronically excited states, is very 
important. This is especially true in the post-discharge phase. Neglecting these aspects 
will lead to non-negligible errors in any model. The results in Fig. [6] warns on the use 
of the Maxwell distribution function for calculating electron-molecule rate coefficients. 
The results reported in Fig. [2l as well as in Ref. [31] can be used only for particular 
conditions i. e. at very high ionization degree when electron-electron Coulomb collisions 
dominates. 

4. Conclusions and Perspectives 

We have presented the dependence of the macroscopic (mean electron energy, vibrational 
temperature) and microscopic (vdf, eedf) quantities in a nitrogen discharge plasma on 
different models for the e-V transitions. As a benchmark model we use the recent 
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database calculated by Laporta et al. Their cross sections are compared with the 
corresponding ones obtained using the analytic expression of Bourdon et al. that 
completes the data-set by scaling-laws the ground state cross sections, and other scaling- 
laws, which restricts the number of interconnecting vibrational levels. Use of these 
restrictive models give results far from the benchmark model while the scaling-law of 
Bourdon et al. gives good agreement, particularly for mean electron energies larger 
than 0.5 eV. Note that, for sake of coherence, all the scaling-laws have been applied 
to the data set of Laporta et al. It should be pointed out that the scaling-law of 
Bourdon et al. has been derived from the accurate set of e-V cross sections of Huo et 
al. This implies that in general complete, ab initio e-V cross sections should be used in 
preference to scaled ones if one wants a good representation of the vibrational kinetics 
of diatomic molecules. In this context we note that complete cross e-V cross section sets 
have recently prepared for CO [56] and O 2 |S7|. In future we plan to use these cross 
sections to extend our approach to the study of more complex mixtures such as air and 
air/hydrocarbon involved in combustion processes, where the presence of excited states 
is very important. 

In these new studies we plan to apply a more realistic form of the applied electric 
held by substituting the present Heaviside-tipe temporal evolution of E/N with shapes 
derived from the experiments [5l|55l[58]- The new results should conhrm the differences 
obtained in the present work in the comparison of models A-C {i.e. on the choice of 
e-V scaling-laws) even though we expect no-negligible variations in the absolute values 
of the relevant macroscopic and microscopic quantities following the selected different 
experimental time dependent voltage forms. 
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